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NOMENCLATURE 

2  , 

a         thermal  diffusivity,   cm  /sec 

a  coefficient  of  linear  thermal  expansion,   1/  K 

A  relative  amplitude  of  electroiiagnetic  radiation 

b  7-coordinate  of  vertex-  of  parabola 

c  x-coordinate  of  vertex  of  parabola 

d  thickness  of  metal  sheet  or  mirror,   cm 

D  diameter  of  mirror,   meters 

f  ratio  of  aperture  to  focal  length 

P  focal  length,   meters 

Fo       Fourier  number,   nondimensional  time 

2        2 

I    intensity,  watts/cm  or  kW/cm 

2      2 

I,    intensity  of  laser  beam,  watts/cm  ,  kW/cm 

k  thermal  conductivity,  watts/cm  K 

I  length  of  objact  expanding 

n  unit  vector  normal  to  surface 

p  distance  from  focus  to  directrix,  meters 

r  reflectivity,  or  radial  coordinate  on  mirror,  meters 

R  characteristic  dimension  of  laser  beam,  or  nondimensional  mirror  radius 

t  temperature,  K 

t  characteristic  temperature,  K 


c 


o. 


t„   initial  temperature,  K 
T    nondimensional  temperature 


9    time,  seconds;  or  angle  between  beam  and  normal,  degrees 

X    coordinate  for  distance  through  sheet  nf  metal  or  mirror,  centimeters, 
meters 

X    nondimensional  distance  through  plate,  x/d 


I.      INTRODUCTION 

Aaong  the  applications  of  high  energy  lasers  is  welding.     For  this  applica- 
tion and  others y  one  would  like  an  estinate  of  the  tiae  required  to  achieve  a 
specified  temperature.     Other  infomation  that  is  useful  is  temperature 
distribution  within  the  metal. 

A  set  of  graphs  has  been  developed  which  pennits  rapid  evaluation  of  time- 
to-temperature,   tempera ttire  profiles,  etc.     First,   the  theory  is  given  for  the 
problem.     Next,  the  limitations  are  outlined  clearly.     Finally,   the  graphs  are 
discussed  and  sample  problems  are  wolfed. 

II.     THEORY  AND  ANAL13I3 

The  time -de pendent,   one-dimensional  equation  for  heat  conduction  is 

8t     8^t  ,tn 

—  "a  — 5  U; 

89     8x^ 
Equation  (1)  is  solved  for  a  plate  of  thickness  d  with  a  constant  heat  fltix  at 
the  surface  x  »  0,  i.e.,  the  front  wall.  The  back  wall,  x  »  d,  is  assumed  to 
be  insulated  so  that 


ax 


(2) 


x-d 


2 

The  constant  heat  flux  is  due  to  a  laser  beam  of  intensity  I,  watts/cm  . 

Consequently,  the  boundary  condition  on  the  front  wall  is 

I  '  -^l^  (3) 


Introduce  the  following  nondiaenBional  yariables  into  equations  (l)  and  (3) 
(a),  X  -  x/d;  (b),  T  =  (t  -  \)/^^i   (c),  t^  -  IdA;  (d)  Fo  •  a  e/d^        U) 


giring 


_8T  ^  afr 

8Fo   8X^ 


(5) 


and 


H  -  -^  <^) 


The  solution  to  the  problea  fonrulated  above  la 
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T-Fo+--(X-~)--%7      ' 5 cos  nrrX  (7) 


n-1 
Dae  to  the  exponential  term  in  the  series,  for  Fo  >  0.^  the  error  in  dropping 
the  series  is  less  than  1  per  cent. 

III.   LIMITATIONS  OF  ANALYTICAL  MDDEL 

The  assumption  of  one-dimensional  heating  is  Talid  so  long  as  R  >  >  d,  >^ere 
R  is  a  characteristic  dimension  of  the  laser  beam.   The  assumption  of  constant 
flux  into  the  metal,  as  stated  by  equation  (3),  may  or  may  not  be  realistic. 
This  point  will  be  discussed  in  more  detail  below. 

Implicit  in  the  equations  (l)  to  (7)  is  the  fact  that  material  properties 
remain  constant.  Thermal  conductivity,  k,  is  a  function  of  temperature.  Like- 
wise thermal  diffusivity,  a  ,  is  a  function  of  temperature. 

Return  again  to  the  question  of  constant  heat  fltuc  due  to  laser  radiation. 
The  laser  beam  has  the  intensity  I, .   The  value  of  intensity  that  is  correct  for 
equation  (3)  is 

I  »  I^jCos  e(l  -  r)  (8) 


Equation  (8)  accounts  for  the  spreading  of  the  radiation  due  to  angle  9 
(See  Fig,  1. )  cmd  the  reflectivity.   The  reflectivity,  r,  ia  a  weak  fiinction 
of  temperature. 

The  charts  are  of  value  only  up  to  the  melting  temperature.  All 
parameters  in  the  problem  vary  rapidly  with  temperature  as  phase  changes  from 
solid  to  liquid.  Furthermore,  the  heat  of  vaporization  is  not  incorporated  in 
the  model.  When  the  surface  becomes  molten,  vapor  forms.   The  vapor  may  be 
absorbing;  as  a  consequence,  the  intensity  of  radiation  may  be  greatly  reduced. 

In  suamary,  the  model  provides  adequate  answers  up  to  the  melting  tempera- 
ture if  the  radiant  heat  flux,  I,  is  correctly  determined  by  equation  (7), 

IV.   DISCUSSION  OF  GRAPHS 

Figures  2  to  4  are  graphs  which  have  been  calculated  using  equation  (7) 
and  the  appropriate  nondimensional  variables,  equations  (4-a)  to  (^d),  which 
were  defined  earlier.   Note  that  there  are  three  graphs,  i, e,,  one  in 
upper  left  hand  (ULH)  comer,  one  in  upper  right  hand  (URH)  corner,  and  one  in 
lower  right  hand  (LRH)  corner.   The  primary  graph  is  in  the  URH  corner;  both 
ULH  and  IRH  graphs  are  for  convenience  in  the  use  of  the  primary  graph.   The 
URH  graph  is  a  plot  of  nondimensional  temperature,  T,  in  terms  of  non- 
dimensional  time,  Fo.   The  Fourier  number  is  essentially  a  nondimensional 
time.  To  use  the  graphs,  one  maat  know  the  Intensity,  I;  the  thickness  of  the 
plate,  d;  the  thermal  conductivity,  k;  and  thermal  diffusivity,  a  ,  A 

characteristic  temperature,  t  ,  is  calculated  using  equation  (4c);  t  is  used 

c  c 

to  enter  the  ULH  graph. 

The  LRH  graph  has  been  drawn  for  aluminum.   It  is  a  plot  of 

(9) 


10 


Figure  1.   Geometry  of  Heating  by  Laser  Radiation, 
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2 
For  aluminum,  the  thermal  diffusivity  is  0.521  cm  /aec.   The  time,  6,  is  used 

as  a  parameter  in  the  LRH  graph. 

Both  the  ULH  and  URH  graphs  are  universal j  i.e.,  they  can  be  applied  to  any 

homogeneous  substance.   However,  the  LRH  graph  is  applicable  only  for  aluminxm. 

A  simple  calculation  can  be  used  to  convert  the  time,  9,  read  from  the  graph 

to  the  time,  0',  for  other  materials.   This  calculation  is  based  on  equation  (9) 

as  follows: 

0'  »  0a  /a' 

The  other  material  has  thermal  diffusivity,  a'.  For  example,  suppose 

2 

0"  «  0.2605  cm  /sec  for  material  "I";  hence 

Q,  .  (1.2  second) (0.521  cm^sec)  ^  3.4  second 
(0.2605  cm  /sec) 

Vlth  the  simple  calculation  outlined  above,  the  graphs  of  Figs.  2,   3,  and  U 

are  universal. 

Some  sample  problems  will  now  be  discussed. 

A.   TEMPERATURE  AT  A  GIVEN  TIME  AND  LOCATION 

Sample  Problem  1 

Assume  a  plate  thickness  of  d  «  1.3  cm  and  an  initial  temperature  of 
tQ  X  300°K.   Assume  I  =  1000  watt/cm^  and  k  =  1. 3  watt/cm°K.   If  time 
9  '  0.4  second  has  elapsed  since  the  laser  beam  was  switched  on,  what  is  the 
temperature  at  front  face? 


Calculate  t  . 
c 


c     (1.3; 
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Enter  the  LRH  graph  with  d  «  1.3  cm.   Move  to  the  right  to  the  curve  for 

6  =  0,U»      Then  move  upward  to  the  curve  for  1=0,  which  is  the  front  face. 

Note  that  Fo  =  0.12  and  T  =  .39.   Move  toward  the  left  imtil  the  line  for 

t  =  1000  K  is  intersected.   Then  move  down  to  the  axis  and  read  t  -  t^  =  390. 
c  0 

Since  the  initial  temperature  was  300  K,  the  front  face  has  a  temperature 
t  =  690°K. 

Sample  Problem  2 

Sample  Problem  2  is  the  same  as  Problem  1  except  that  the  temperature  at 

X  =  0.52  cm  is  desired.  Calculate  X  =  x/d  =  0.52/1.3  «  0.4.  Where  the 

vertical  line  running  between  the  LRH  and  URH  graphs  intersects  I  =  0,U   curve, 

draw  a  horiaontal  line.   At  the  intersection  with  the  t  =  1000  K  curve, 

c  ' 

descend  vertically  to  the  axis.  Read  t  -  t-  «  115»  At  a  distance  x  =  0.52  cm 
and  a  time  9  »  0.4  sec,  the  temperature  is  4-15  K. 

B.   TIME  TO  MELT  FRONT  SURFACE 

Sample  Problem  3 

The  melting  temperature  of  aluminum  is  933  K.   Assume  t-  »  333  K. 

If  a  beeun  has  an  intensity  of  2000  watt/cm  ,  how  much  time  is  required  to 

begin  melting  the  front  surface  of  an  altuninum  plate  with  thickness  0.65  cm? 

CaJculate  t  . 
c 

t  .  (20O0K0,6?)  .     Oj 

C         X,  ) 

Subtracting  333  K  from  933  K,  one  enters  the  ULH  graph  at  600  K  and  moves 
upward  to  t  =  4000  K  curve.   Moving  horizontally,  one  crosses  the  T-axis 
at  0.152  and  stops  at  the  intersection  of  I  =  0  curve.   Then  a  line  is  dropped 
vertically  crossing  the  Fourier  number  axis  at  0.018.   Draw  a  horizontal  ciirve 
in  the  LRH  graph  at  d  «  0,65  cm.   The  intersection  of  the  horizontal  and  vertical 
lines  shows  that  0  =  .015  second.  Note  that  the  front  face  is  ready  to  melt  while 
the  ten^jerature  of  the  back  face  (X  «  l)  remains  t  «  333°K. 
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V.     MIMtOR  DISTORTION  DTffi  TO  RADIANT  HEATING 

High  powered  lasers  often  use  metal  mirrors  which  undergo  thermal  expansion 
and  optical  distortion  due  to  absorption  of  radiation.      Depending  on  mirror 
geometry,    complex,   three-dimensional  calculations  may  be  required.      In  this  section, 
a  sample  calculation  for  determining  one -dimensional  mirror  distortion  is  illustrated. 

The  fornula  for  themal  one-diiaensional  strain  is 

^  .  AAt  (10) 

where  a  is  the  coefficient  of  linear  thermal  expansion,  t   is  the  length  of 

the  object  expanding,  A-t  is  the  increase  in  its  length,  and  At  is  the  temperature 

change  involved. 

The  linear  thernal  expansion  coefficient,  a,  is  a  function  of  temperature. 
Table  I  gives  Its  values  for  different  temperatures  of  aluainua.  Values  for 
other  netals  as  a  function  of  temperature  can  be  found  In  the  American  Instituts 
of  Physics  Handbook  ^f .  l/. 

Table  I,   Temperature  Variation  of  Linear  Thermal 

Expansion  Coefficient  for  Aluminum  /fief.   l/. 

t  (°K) a  (10"^ A) 

300  23.2 

350  2ii.l 

^00  2i^.9 

500  26.4 

600  28.3 

700  30.7 
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The  expansion  is  given  by  equation  (11). 


a 
^         J 


a(t)(t  -  tQ)dx  (11) 


0 

where  x   is  the  coordinate  for  the  distance  through  the  mirror  and  d  is  the 
mirror  thickness.  For  this  sample  problem,  numerical  evaluation  has  been 

chosen.  The  following  characteristics  for  the  mirror  were  used:  parabolic 

2 
mirror,  mirror  radius  «  0.05  m,  f  number  =  1.0,  intensity  of  beam  =  4-0  kW  cm  , 

mirror  absorption  »  0.01,  and  Fresnel  number  =  10, 

Using  the  facts  that  the  mirror  is  jaarabolic  with  a  f  number  of  1.0  and 

thickness  at  the  center  of  2  cm,  the  thickness  of  the  mirror  at  other  radial 

points  can  be  calculated, 

-  ^  Focal  length  ^  F 
Diameter     D 

F  =  Df  =  (0.1  m)(l)  =  0.1  ■ 

The  standard  form  for  a  parabola  is 

-   2p(y  -  b)  -  (r  -  c)^  (12) 

where  p  is  the  distance  from  the  directrix  to  the  focus  of  the  parabola  and 

(c,b)  are  the  coordinates  of  its  vertex.  Thus  p/2  is  thA  focal  length  of  the 

mirror.  So  p/2  =  0.1  m,  and  p  =  0.2  m.  Putting  the  vertex  (c,b)  of  the 

parabola  at  the  origin  yields  c  «  b  «  0.  Therefore,  the  equation  of  th« 

2  2 

parabola  is  2(0.2)(y  -  O)  «  (r  -  O)  or  y  »  2.5r  .  Solving  this  equation  for 

different  values  of  r  gives  the  increase  in  the  thickness  of  the  mirror  at 

different  radii.  The  actual  thickness  can  be  found  by   adding  2  cm,  which  is  the 

thickness  at  the  center  of  the  mirror.  The  thickness  of  the  mirror  at  different 

radial  locations  is  given  in  Table  II, 
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Table  II.   Mirror  Thickness,  Intensity  of  Absorbed  Power, 
and  Characteristic  Temperature 


R 

r,  m 

d,  cm 

Intensity, 

W/cm^ 

v'^ 

0 

0.0 

2.000 

AOO 

615 

0.3 

0.015 

2.056 

296 

468 

0.7 

0.035 

2.306 

77 

137 

1.0 

0.05 

2.625 

5 

9.77 

To  find  the  intensity  of  the  laser  radiation  at  rarioua  radii  on  the  mirror, 

the  amplitude  profile  from  Figure  3.11  from  Lengyel  ^ef.  2j  was  used.   The  absorbed 

2 
intensity  at  the  center  of  the  mirror  was  assumed  to  be  4-00  watts/cm  which 

was  obtained  by  multiplying  radiation  intensity  \jj  mirror  absorption 

(40kW/c»  X  O.Ol).   Amplitude  and  intensity  are  related  by 


I, 


^0 


(13) 


'0 

where  Aq  and  A-,  are  the  amplitudes  at  two  different  points  (from  the  referenced 

graph),  and  I-  is  the  known  intensity  at  the  mirror  axis.   The  intensities  also 

are  tabulated  in  Tabl«  II. 

In  this  problem,  one-dimensional  heating  is  assumed;  i.e.,  the  effects  of 

heat  transfer  in  the  radial  direction  are  ignored.  To  solve  the  problem,  it  is 

necessary  to  calculate  t^.   This  can  be  done  using  equation  (4c)  where 

k  =  1.3  watt/cm  K.   The  values  of  t  for  different  values  of  d  are  listed  in 

c 

Table  II. 

19 


To  evaluate  the  Integral  numerically,  the  mirror  is  imagined  as  being 
divided  into  slabs.   The  expansion  of  each  of  these  slabs  is  calculated 
based  on  the  mean  slab  temperature.   These  calculations  can  be  made  using 
the  set  of  graphs,  Figs.  2  to  4,  explained  earlier  in  this  paper.  The 
calculation  methods  become  apparent  by  referring  to  Table  III  which  svuranarizes 
the  results  for  the  centerline  position  on  the  mirror. 

The  graph.  Fig.  5,  is  a  profile  of  the  thermal  expansion.  From  the 
graph,  it  can  be  seen  that  considerable  distortion  will  occur,  even  at  times 
as  small  as  0.3  sec.   Modest  quality  optics  have  a  tolerance  less  than  one 
quarter  of  the  wavelength  of  the  light.   The  distortion  in  this  problem 
would  be  considerably  more  than  a  quarter  wavelength. 

VI.   CONCLUSIONS 

Graphs  have  been  developed  for  the  solution  of  many  types  of  problems 
involved  with  laser  heating.  Rapid  solutions  become  possible  by  suitable 
use  of  the  figures.  Sample  problems  were  presented  for  temperature  profiles 
within  the  plate,  time  to  melt  the  front  surface,  and  mirror  distortion  due  to 
thermal  expansion. 
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Table  III.     Svunmary  of  Calculations  for  Thermal  Expansion  at  Mirror  Center 


X  =  x/d  X,   cm  t  -  t^*  <t  -  t^>  a  Ax  Sx,   microns 


0.6  1.2 


UO 


0.1  0.2  80 


0.2  0.4  50 


0.3  0.6  25 


O.U  0.8  10 


0.8  1.6  0 


110  .252       .2  5.5U 


65  .2U       .2  3.172 


37.5  .239       .2  1.793 


17.5  .235       .2  0.823 


6  .233       .A  0.559 


.232       .4.  0.093 


1.0  2.0  0 

Total  Thermal  Expansion,  At  11.983  microns 

^Obtained  from  Fig.    3  with  d  =  2  cm  and  t     =  615°K. 

c 

<t  -  t-.>  average  slab  temperature  in  excess  of  initial  temperature  t.,  K. 

a  coefficient  of  linear  thermal  expansion,  microns/cm  K. 

Ax  slab  thickness,  cm. 

8x  thermal  expansion  of  slab,  microns. 
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Figure  5.   Profile  of  Thermal  Expansion  at  300  Mllliseca, 
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